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ABSTRACT 

While the star formation rates and morphologies of galaxies have long been known to correlate 
with their local environment, the process by which these correlations are generated is not well un- 
derstood. Galaxy groups are thought to play an important role in shaping the physical properties 
of galaxies before entering massive clusters at low redshift, and transformations of satellite galaxies 
likely dominate the buildup of local environmental correlations. To illuminate the physical processes 
that shape galaxy evolution in dense environments, we study a sample of 116 X-ray selected galaxy 
groups at z = 0.2 — 1 with halo masses of 10^^ — 10^^ M© and centroids determined with weak lensing. 
We analyze morphologies based on HST imaging and colors determined from 31 photometric bands 
for a stellar mass-limited population of 923 satellite galaxies and a comparison sample of 16644 field 
galaxies. Controlling for variations in stellar mass across environments, we find significant trends 
in the colors and morphologies of satellite galaxies with group-centric distance and across cosmic 
time. Specifically at low stellar mass (log(M^/M0) = 9.8 — 10.3), the fraction of disk-dominated star- 
forming galaxies declines from > 50% among field galaxies to < 20% among satellites near the centers 
of groups. This decline is accompanied by a rise in quenched galaxies with intermediate bulge+disk 
morphologies, and only a weak increase in red bulge-dominated systems. These results show that 
both color and morphology are influenced by a galaxy's location within a group halo. We suggest that 
strangulation and disk fading alone are insufficient to explain the observed morphological dependence 
on environment, and that galaxy mergers or close tidal encounters must play a role in building up the 
population of quenched galaxies with bulges seen in dense environments at low redshift. 
Subject headings: galaxies: bulges - galaxies: clusters: general - galaxies: evolution - galaxies: halos 
- galaxies: statistics - X-rays: galaxies: clusters 



1. INTRODUCTION 

Galaxy properties are correlated with their environ- 
ment. Dense regions host galaxies with greater masses, 
lower star-formation rates, and morphologies that are 
more bulge-dominated than in low density regions. Since 
the early work of Dressier (1980), studies of environmen- 
tal correlations have expanded to a variety of indicators 
of star formation, morphology, and environment. See 
Blanton & Moustakas (2009) for a recent review of these 
correlations in the local universe. 

Numerous processes may be culpable for the depen- 
dence of galaxy properties on environment (see Boselli & 
Gavazzi 2006 for a review). Galaxy interactions through 
mergers or tides can disrupt stellar kinematics and re- 
move gas. These interactions can also apply torques that 
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drive gas inward, perhaps feeding star formation or a 
central black hole. Halos can play a role through tidal 
forces and dynamical friction. Pressure from hot dense 
gas inside halos may strip gas from inf ailing galaxies, and 
shock heating in massive halos can prevent accretion of 
cold gas that would feed star formation. 

Despite extensive observational and theoretical work, 
the dominant physical mechanisms responsible for the 
environmental correlations remain unclear. An impor- 
tant clue is that the scale on which these correlations ap- 
pear is similar to the size of the dark matter halos hosting 
galaxies (e.g., Kauffmann et al. 2004; Blanton & Berlind 
2007). Relatedly, satellite galaxies have a primary role 
in the buildup of the red sequence in dense environments 
(e.g., Weinmann et al. 2006; van den Bosch et al. 2008; 
Wetzel et al. 2012a; Peng et al. 2012). 

Since many galaxy characteristics are interrelated, it 
is important to study the correlations for each prop- 
erty independently while fixing other variables. When 
constraining environmental effects, one ideally controls 
for differences in redshift, stellar mass, halo mass, and 
location within a group. The colors and morphologies 
of galaxies are also correlated, so it is advantageous 
to split galaxies simultaneously by color and morpho- 
logical classes to distinguish between processes that af- 
fect star formation rates and structural properties differ- 
ently. Several studies have suggested that environmental 
processes affect star formation more significantly than 
morphology (e.g., Kauffmann et al. 2004; Blanton et al. 
2005; Christlein & Zabludoff 2005; Weinmann et al. 2009; 
Kovac et al. 2010), suggesting that the gas that feeds star 
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formation is depleted or removed without significantly al- 
tering galactic structure or stellar kinematics. 

While many of the observational studies have focused 
on large surveys at low redshift or on massive clusters 
at higher redshifts, less is known about the more com- 
mon group-scale environments and their evolution over 
time. In this paper, we study the colors and morpholo- 
gies of galaxies in groups spanning the redshift range 
z = 0.2 — 1. We focus in particular on the properties 
of satellite galaxies as a function of stellar mass, group- 
centric distance, and redshift. This sample of groups 
has been identified from the COSMOS field (Scoville 
et al. 2007) based on their extended X-ray emission and 
have masses determined with weak lensing in the range 
10^^ - 10^"^ Mq (Leauthaud et al. 2010). Galaxy mem- 
bership has been assigned using precise photometric red- 
shifts with a Bayesian procedure tested extensively with 
mock catalogs and a spectroscopic subsample (George 
et al. 2011). Furthermore, by finding galaxies that maxi- 
mize the lensing signal on small scales we have identified 
central galaxies that accurately trace the center of mass 
of their halos (George et al. 2012). Additional studies 
of these groups with magnification lensing (Ford et al. 
2012; Schmidt et al. 2012) and clustering (Allevato et al. 
2012) have confirmed the halo mass estimates. These in- 
gredients make for a robust characterization of group en- 
vironments with which to study satellite properties and 
their evolution over time. 

We describe the data used in this analysis in Section 2 
along with mock catalogs used to estimate the purity 
of the satellite population analyzed. Section 3 presents 
the results including our primary finding of a declining 
fraction of blue late disk galaxies and a rise in red early 
types toward group centers and over cosmic time. In 
Section 4, we discuss the implications of these findings 
for the physical mechanisms that shape galaxy evolution, 
and put this work in the context of previous studies and 
possibilities for the future. 

2. DATA AND MOCKS 

The COSMOS field (Scoville et al. 2007) has been the 
subject of a broad array of mult i- wavelength imaging 
and spectroscopy campaigns. In this study we make use 
of several pieces of this data: galaxy colors, photomet- 
ric redshifts, and stellar masses derived from over thirty 
bands of ultraviolet, optical, and infrared data (Ilbert 
et al. 2009, 2010; Bundy et al. 2010); galaxy morphologies 
determined by Scarlata et al. (2007) from high-resolution 
imaging taken by the Advanced Camera for Surveys 
(ACS) aboard the Hubble Space Telescope (Koekemoer 
et al. 2007); and a catalog of massive galaxy groups se- 
lected based on their X-ray emission seen with XMM- 
Newton and Chandra and for which halo mass profiles 
have been characterized using weak lensing with Hub- 
ble imaging (Leauthaud et al. 2010; George et al. 2011, 
2012). Most of these data have been compiled in the 
group member catalog published in George et al. (2011) 
and we refer the reader to that paper and references 
therein for further details including the flux and mass 
limits for galaxies and groups. Here we briefly explain 
the stellar masses, colors, and morphologies used in this 
paper, as well as the identiflcation of group centers and 
the mock catalogs used for validating group membership 
assignment and estimating the effects of contamination 



from fleld galaxies. 

2.1. Group Membership and Mocks 

As a measure of galaxy environment, we use a cata- 
log of groups described in George et al. (2011). These 
groups have been identifled through a wavelet detection 
of extended X-ray emission following Finoguenov et al. 
(2010) and the sample is an expanded version of an earlier 
COSMOS catalog (Finoguenov et al. 2007) with deeper 
X-ray data. A Bayesian membership algorithm is used 
to associate galaxies with groups based on precise pho- 
tometric redshifts and proximity to the X-ray position, 
with a prior accounting for the relative fleld population as 
a function of magnitude and redshift. Following George 
et al. (2011), we select galaxies with membership prob- 
ability Pmem > 0.5 amoug the clean sample of groups, 
avoiding masked regions, groups with fewer than four 
members, overlapping groups, and possible projections 
(flag_include=1 in the catalog). Halo masses for these 
groups have been measured with stacked weak lensing to 
be in the range M200C ^ 10^^ — IO^^Mq and are correlated 
with X-ray luminosity (Leauthaud et al. 2010; George 
et al. 2012). Over this limited range in halo masses, we 
have not detected signiflcant trends in galaxy properties 
with X-ray luminosity, but to avoid the influence of out- 
liers we eliminate a small number of groups with halo 
masses estimated from their X-ray luminosities follow- 
ing Leauthaud et al. (2010) that are above or below the 
10^^ - 10^"^ Mo range. 

We split galaxies into central, satellite, and fleld pop- 
ulations across several stellar mass and redshift bins to 
isolate dependences among different properties. For each 
group, the central galaxy is deflned as the most mas- 
sive group member within a projected distance of the 
X-ray position equal to the scale radius of a Navarro- 
Frenk- White (Navarro et al. 1996) proflle. George et al. 
(2012) showed with weak lensing tests that this deflnition 
traces the halo center to within roughly 75 kpc, though 
in roughly 30% of groups the center is still ambiguous. In 
this paper we do not speciflcally study the central galax- 
ies in these X-ray selected halos because their abundance 
limits statistical constraints on the properties of this pop- 
ulation. The remaining group members are called satel- 
lites, and are the primary focus of this work. Field galax- 
ies are taken to be those not assigned to any extended 
X-ray source and are thus expected to reside in halos less 
massive than the X-ray detection limit (Pmem = from 
George et al. 2011; see Figure 1 of that paper for the X- 
ray flux limit as a function of redshift). Because we do 
not detect faint X-ray groups at high redshift, the halo 
mass range occupied by the fleld population can evolve 
due to this selection effect. However, at low redshift only 
~ 10% of galaxies in our sample live in X-ray groups, so 
we expect the evolving halo mass limit to have only a 
small effect on the statistics of the fleld population. Ta- 
ble 1 gives the size of each of these samples used in our 
analysis. 

Extensive tests of the accuracy of the photometric red- 
shifts and membership assignment algorithm have been 
carried out using mock catalogs from simulations as well 
as a subsample of galaxies with spectroscopic redshifts 
(George et al. 2011). The overall purity and completeness 
of the member selection within 0.5i?200c is 84% and 92%, 
respectively, down to our flux limit of F814W = 24.2. 
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TABLE 1 
Environmental Census 



Stellar Mass [log(M^/M0)] 
Type [9.8,10.3) [10.3,10.8) 



z = 0.2- 


0.5; A^groups = 47 




Satellites 


237 


218 


Field 


2455 


1993 


z = 0.5- 


0.8;Argroups = 40 




Satellites 


122 


189 


Field 


4497 


3529 


z = 0.8- 


1.0;A^groups = 29 




Satellites 




157 


Field 




4170 


Note. 


— No contamination 


correc- 


tions have been applied. Ellipses 


denote 



a bin below our stellar mass completeness 
limit. 

The accuracy of member selection depends most signifi- 
cantly on distance from the group center and on the flux 
of a galaxy. The first effect is due to the decreasing den- 
sity of true members with group-centric distance, which 
implies that selecting galaxies uniformly with radius (out 
to i^200c) will result in higher contamination from field 
galaxies in the outskirts. The second effect from galaxy 
flux arises from the decreasing precision of photomet- 
ric redshifts for fainter galaxies, ranging from ^ 0.01 
at F814W < 22.5 to a, = 0.03 at F814W = 24. It is 
notable that with the many filters used in constructing 
these photometric redshifts, precision does not depend 
significantly on spectral type, so red and blue galaxies 
of a given magnitude are selected with similar complete- 
ness. 

In our analyses of member populations, we correct for 
contamination from field galaxies using the mock cata- 
logs described in George et al. (2011). A halo catalog is 
constructed from a LasDamas simulation (C. K. McBride 
et al., in preparation)^ and populated with group mem- 
bers following the halo occupation model of Leauthaud 
et al. (2012). Mock galaxies are matched to real COS- 
MOS galaxies in narrow bins of stellar mass and redshift 
to assign magnitudes in the F814W band, which in turn 
are used to assign mock photometric redshift errors. We 
estimate the purity of the member selection algorithm 
in bins of group-centric distance and magnitude, and as- 
sume that the colors and morphologies of contaminating 
field galaxies are representative of the field population at 
that magnitude. Contamination corrections are applied 
to the satellite populations only, since the field popula- 
tion is much larger. These corrections only account for 
contamination of the member list and not for incomplete- 
ness; the latter effect is smaller, does not significantly 
depend on group-centric distance, and should be fairly 
uniform across galaxy types since the precision of photo- 
metric redshifts is effectively achromatic (see Figures 2 
and 5 of George et al. 2011). 

2.2. Stellar Masses, Colors, and Morphologies 

^ Details regarding the LasDamas "Consuelo" simulation 
can be found at http://lss.phy.vanderbilt.edu/lasdamas/ 
simulations . html 



Stellar masses and spectral classes have been deter- 
mined by fitting stellar population synthesis models to 
the spectral energy distributions of galaxies, varying the 
age, amount of dust extinction, and metallicity in the 
models. We use the stellar masses from Bundy et al. 
(2010) and color classifications from Ilbert et al. (2010), 
both of which are based on Bruzual & Chariot (2003) 
models fit with a Chabrier (2003) initial mass function 
but take somewhat different approaches to fitting the 
data. Using the unextincted rest-frame color (NUV — r+) 
of the best-fitting template for each galaxy, we split them 
into two color classes: 

• red: (NUV - r+) > 3.5 

• blue: (NUV - r+) < 3.5. 

Ilbert et al. (2010) showed that this color cut corresponds 
with a specific star formation rate of roughly 10~^^yr~^. 
It is worth emphasizing that these colors are corrected 
for dust extinction, so that the red sample selects a pop- 
ulation of quiescent galaxies similar to those identified 
using two color cuts based on near-ultraviolet, optical, 
and near-infrared measurements, minimizing contamina- 
tion from dusty star- forming galaxies (e.g., Bundy et al. 
2010). 

Galaxies are classified morphologically based on au- 
tomated measurements of several structural parameters 
from ACS imaging in the F814W band. We use the 
Zurich Estimator of Structural Types (ZEST; Scarlata 
et al. 2007) catalog, which is derived from measurements 
of the asymmetry, concentration, Gini coefficient, second 
moment of the brightest pixels, and ellipticity. Scarlata 
et al. (2007) performed a principal component analysis of 
these measured quantities to identify the most important 
structural parameters, and divided galaxies into classes 
(early type, disk, and irregular) based on their location 
in this principal component space. A further subdivi- 
sion of the disk category was derived by correlating mea- 
surements of the Sersic index for a bright sample with 
Iab < 22.5 with location in the principal component 
space. This parametrization provides a rough indication 
of the bulge-to-disk ratio, with four classes ranging from 
bulge-dominated to bulgeless disks. For statistical pur- 
poses, we consolidate some of these classifications into 
three morphological categories: 

• early type spheroidals and bulge-dominated disks 
including even relatively inclined SO galaxies 
(ZEST types 1 and 2.0) 

• intermediate bulge+disk galaxies (ZEST type 2.1) 

• late disks with little or no bulge component (ZEST 
types 2.2 and 2.3). 

Irregulars and unclassified galaxies make up a small frac- 
tion of the sample, typically not more than a few percent 
in any stellar mass bin. We note that morphological 
K-corrections, which have not been applied, should not 
have a large effect on our sample because the F814W 
band probes rest-frame optical wavelengths over our en- 
tire redshift range (e.g., Lotz et al. 2004; Cassata et al. 
2005; Bundy et al. 2010). 
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Fig. 1. — Group members as a function of stellar mass and distance from the group center. Colors for each galaxy represent the average 
unextincted rest-frame template (NUV — r+) color, with shading proportional to the logarithmic surface brightness fi. The gray band at the 
bottom of the high-redshift panels shows the stellar mass completeness limit for a passive population calculated with our flux limits F814W 
= 24.2 and Kg = 24 following the approach of Bundy et al. (2010) (see also Figure 1 of George et al. 2011). Central galaxies in these groups 
are shown in the light gray band on the left side of each of the outer panels. The middle frame shows morphological classiflcations for a 
random sample of galaxies chosen to span the range of colors observed; objects in this panel are sorted vertically by color and horizontal 
offsets within each classiflcation are arbitrary. 
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3. RESULTS 

We begin our analysis with a visualization of a few 
properties of the galaxies in our sample. Figure 1 shows 
the distributions of stellar masses and group-centric dis- 
tances in different redshift bins for all galaxies selected as 
group members. Each point represents one object and is 
displayed using the ACS image of the galaxy (Koekemoer 
et al. 2007), colored according to the average unextincted 
rest-frame template (NUV — r+) color. Centrals are po- 
sitioned on the left of each panel with small horizontal 
offsets, and satellites are plotted at their distance from 
the corresponding central. Each image is basically the 
set of adjacent pixels with flux above a noise threshold 
shown with a logarithmic surface brightness scale that 
has its maximum set to the peak value for each galaxy. 
We note that the apparent size of a galaxy on the plot is 
quite sensitive to its flux since images of brighter galax- 
ies have more pixels above the noise threshold, and this 
is not a perfect indicator of the physical effective radius 
of a galaxy. A small fraction of objects have deblend- 
ing issues or edge effects from the cutout image size, but 
we note that the images are processed independently for 
visualization and analysis purposes. 

Several trends are evident when visualizing galaxies in 
this manner. Stellar mass is a strong determinant of 
galaxy properties; massive galaxies are more likely to be 
large, red, and spheroidal. Physical properties also de- 
pend on the location within a group. Central galaxies 
are massive (by definition), but also typically red and 
elliptical. Blue centrals tend to be less massive than red 
centrals at high redshifts even within the fairly narrow 
range of halo masses studied here, as discovered by Tin- 
ker et al. (2012). Satellites closer to group centers are 
more likely to be red and show fewer spiral features than 
in the outskirts, particularly at low stellar masses. There 
is a relative dearth of low mass satellites near group cen- 
ters; this could be a hint of satellite depletion due to 
mergers, but challenges in measuring photometry of faint 
objects near massive extended central galaxies could be 
a contributing factor so further investigation is needed. 
Similar evidence of mass segregation and the influence of 
mass and environment on star formation has been seen 
in groups and clusters out to z ^ 1 (Muzzin et al. 2012; 
Presotto et al. 2012). Many of these trends can be seen 
across the entire redshift range studied here, indicating 
that both stellar mass and halo environment play a role 
in determining galaxy properties at least since z = 1. 
Sample variance due to the finite survey volume does af- 
fect our ability to measure absolute redshift trends as 
the number densities vary significantly due to large-scale 
structures (particularly evident at z = 0.3 — 0.4), but the 
relative fractions of different populations should be less 
affected. 

The central panel of Figure 1 shows a few examples of 
the three morphological classifications from ZEST, cho- 
sen randomly to span the range of colors seen from red- 
shifts z = 0A- 0.6 and with log(M^/M0) > 10. While 
these classifications do not always agree with one's vi- 
sual impression, there are clear differences in structural 
parameters among the classes, and correlations between 
these automated measurements and a galaxy's position 
in a group can provide an interesting test of the depen- 
dence of morphology on environment. We emphasize 



that the ZEST morphologies are correlated with, but not 
identical to, traditional visual classification of ellipticals 
and spirals. We compare results from multiple morpho- 
logical indicators in Section 3.3. 

3.1. Radial Trends: Blue Late Disks into Red 
Bulge-h disks 

We can quantify some of the trends from Figure 1 by 
measuring the fraction of galaxies of a given color and 
morphology as a function of stellar mass, group-centric 
distance, and redshift. Figure 2 shows the fraction of 
galaxies in each of the six combinations of color and mor- 
phology categories described in Section 2.2. For exam- 
ple, the cyan triangles represent the fraction of galaxies 
in the stellar mass and redshift range shown that are 
both blue and have late disk morphology. This popula- 
tion makes up the majority among field galaxies (shown 
at > i^200c) but its fraction declines among satel- 
lites, contributing less than 20% of the satellite popu- 
lation dX R < i^200c/3. Meanwhile, the proportion of 
red bulge+disk galaxies rises from 7% in the field to 40% 
among satellites in the inner radial bin. 
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Fig. 2. — Color and morphological fractions as a function of 
group-centric distance. These fractions are calculated for satellites 
in three equally-spaced bins out to i^200c after applying contam- 
ination corrections. The field population is plotted to the right. 
Points are assigned small horizontal offsets for clarity. Error bars 
are the la standard deviation of 500 bootstrap samples. There is a 
clear transition from blue late disks dominating among field galax- 
ies and outer satellites to red bulge+disks among inner satellites. 

Figure 2 highlights the most significant environmental 
trends in our sample by focusing on low mass, low red- 
shift galaxies. We repeat the exercise in Figure 3 with 
higher mass and redshift bins to study how these trends 
vary. The broad picture is similar; blue late disks dom- 
inate in the field while the red bulge+disk population 
becomes more prominent toward group centers. Among 
low mass galaxies at z > 0.5, the blue late disks dom- 
inate everywhere, while more massive galaxies at lower 
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redshift have a substantial population of red spheroidals. 
Red late disks and blue spheroidals do not make up a 
large portion of galaxies at any mass or redshift studied. 
The red fraction can be determined from these plots by 
summing the three red lines, and similarly the spheroidal 
fraction is the sum of the solid lines with circular mark- 
ers. The red fraction rises toward group centers for low 
mass galaxies, but is relatively flat among massive galax- 
ies. We note that the abscissa for these plots is the pro- 
jected group-centric distance, and that the true radial 
trends measured in spherical shells are likely more sig- 
nificant than observed. 
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Fig. 3. — Color and morphological fractions of group members as 
a function of distance from the group center, for different bins of 
stellar mass (columns) and redshift (rows). Line styles and error 
bars are as defined in Figure 2. Figure 2 is repeated in the top 
left panel, for comparison with weakening environmental trends at 
higher mass and redshift. 



3.2. Redshift Trends 

Since satellites tend to fall toward halo centers, group- 
centric distance is related to the timescale that galaxies 
have been inside the group. The range of redshift s sam- 
pled with this data set provides another measure of time 
to study evolution. Figure 4 is a transpose of Figure 3 to 
show the redshift trends in the distribution of colors and 
morphologies. Again, there is a decline in blue late disks 
among low-mass satellites near the centers of groups, now 
compensated by a rise in both red spheroidals and red 
bulge+disks. This trend weakens away from group cen- 
ters and at higher masses. 

3.3. Checks for Systematics 
3.3.1. Environment 
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Fig. 4. — Color and morphological fractions as a function of red- 
shift, for different bins of stellar mass (columns) and environment 
(rows). Inner and outer satellites are separated at a projected 
group-centric distance of 0.5i?200c- Line styles and error bars are 
as defined in Figure 2. 

There are several possible biases or other effects in the 
data to consider in order to ensure the robustness of these 
results. First, we revisit the contamination of our satel- 
lite sample due to interloping field galaxies, discussed 
in Section 2.1. In Figures 2, 3, and 4 we have plotted 
values of population fractions for satellites corrected for 
contamination estimated from mock catalogs. Contami- 
nation corrections are always smaller than the statistical 
error bars estimated via bootstrapping, except for low 
mass blue disk galaxies where it is 10% greater than the 
error because the field population is so large. Though the 
sample of satellites near i?200c is significantly contam- 
inated by field galaxies, the corrections to the relative 
fractions of each galaxy type are small because the field 
populations are not markedly different from the outer 
satellites. 

3.3.2. Color Distribution 

When classifying galaxies by their spectral energy dis- 
tributions (SEDs), our primary aim is to distinguish star- 
forming galaxies from those that are quenched. Though 
traditional indicators from emission lines or spectral 
breaks are not directly available from photometry, the 
31-band SEDs used here provide a wealth of information 
about spectral types, including an estimate of dust ex- 
tinction that separates star-forming galaxies that appear 
red due to dust from those that are truly passive. The 
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Fig. 5. — Distribution of rest-frame, extinction-corrected, tem- 
plate colors by morphological type. Galaxies from all environments 
are included. Vertical dotted lines show the cut used to segregate 
red and blue galaxies. 

template-based, extinction-corrected (NUV — r+) colors 
used in this paper generally have a bimodal distribution 
with the color cut from Section 2.2 falling on the red end 
of the "green valley." The color distributions for differ- 
ent morphological types are shown in Figure 5. Shifting 
the cut slightly in either direction shifts the amplitude 
of the red fraction in Figures 2, 3, and 4 up or down, 
but the trends with group-centric distance and redshift 
do not vary significantly. We have tested an alternative 
"red sequence" selection suggested by Ilbert et al. (2010), 
NUV - r+ > O.51og(M^/M0) - 0.8z - 0.5, using rest- 
frame absolute magnitudes but applying no extinction 
correction. The transformation of blue late disks into 
red bulge+disks among low mass satellites is still evident. 
Similarly, applying the two-color cut in the NUV — r+ , 
r+ — J plane used by Bundy et al. (2010) does not qual- 
itatively change our results. 

3.3.3. Morphologies 

Morphological classification is a challenging problem 
and significant scatter exists between the types assigned 
to galaxies in both visual and automated analyses. In 
general, visual analyses emphasize the presence or ab- 
sence of spiral features, categorizing objects as spirals 
or ellipticals, often with an intermediate class of SOs 
grouped with ellipticals. Automated analyses measure 
structural parameters such as concentration and asym- 
metry which are then generally tied to a training set of 
visual classifications. The correlation between the prop- 
erties measured is imperfect, so we test the impact on 
our results of using a variety of automated morpholog- 
ical classifications, all based on the ACS F814W imag- 
ing. The alternative classifications come from Tasca et al. 
(2009), which presented three separate techniques. 

The differences between the results of each catalog and 
those from the ZEST classification used in this paper are 



driven by how bulge+disk galaxies are classified, since 
the other catalogs do not split the spiral/disk category 
into multiple bins as ZEST does. For instance, among 
the 237 satellites used for Figure 2, the ZEST classifi- 
cations are 22% spheroidal, 36% bulge+disk, 37% late 
disk, and 5% irregular or unclassified. The three sep- 
arate classifications for these galaxies from Tasca et al. 
(2009) vary between 39 - 59% E/SOs and 59 - 40% spi- 
rals. When bulge+disk galaxies are mostly classified as 
E/SOs, the spheroidal populations (both blue and red) in 
Figures 3 and 4 are significantly elevated but show sim- 
ilar trends with group-centric distance and redshift. In 
the opposite case where bulge+disk galaxies are mostly 
treated as disks, the spheroidal fraction is essentially un- 
changed. While the dominant population in each panel 
of Figure 3 and 4 can change depending on which cate- 
gory the intermediate bulge+disk galaxies fall into, the 
radial trend in low mass satellites is unchanged; the blue 
late type fraction declines toward group centers and is 
compensated with a rise in red early types. The redshift 
trend for this transition is strongest for low mass satel- 
lites near group centers and weakens toward larger radii 
and stellar masses. Though the scatter between mor- 
phological classifications signals that our results should 
be interpreted and compared to others with caution, the 
significant trends with group-centric distance and red- 
shift for the intermediate bulge+disk galaxies suggests 
that the ZEST classification has identified a population 
in transition. 

3.3.4. A Population of Blue Spheroidals 

While the tests described above suggest that our mea- 
sures of environment, color, and morphology are ro- 
bust, we do note a puzzling population of massive blue 
spheroidals that can be seen most clearly in the bot- 
tom right panel of Figure 4 as well as the right col- 
umn of Figure 5. Those plots suggests that half of all 
spheroidals in that stellar mass range are blue, exceed- 
ing measurements in other studies (e.g. Kaviraj et al. 
2007, 2008; Bamford et al. 2009; Schawinski et al. 2009; 
Ilbert et al. 2010), but see also Cross et al. (2004) who 
found a large fraction of blue ellipticals in a luminosity- 
selected sample at moderate redshift. We see a similarly 
large blue fraction among spheroidals at higher masses 
(log(M^/M0) > 10.8, not plotted) where spheroidals 
make up a higher proportion of all galaxies, and most 
prominently at low redshift {z < 0.5). Visual inspection 
of the ACS images of these galaxies suggest that 30—40% 
show spiral or irregular features, although the structural 
parameters measured by the automated classifiers are all 
consistent with the red spheroidal population. We have 
also investigated the UV-optical and optical-IR colors of 
these massive blue spheroidals prior to template fitting, 
in addition to optical spectra for 30 of these objects de- 
scribed in George et al. (2011). These data suggest that 
blue spheroidals are not typical star-forming galaxies but 
may have a small amount of recent star formation to 
which rest-frame UV measurements are particularly sen- 
sitive (e.g. Kaviraj et al. 2007, 2008). Another possible 
explanation is that there are unusual stellar populations 
that are not represented in the template fits used to de- 
rive our (NUV — r+) colors (see e.g. Smith et al. 2012, 
for a discussion of galaxy properties contributing to the 
"UV upturn"). While this population is interesting in 
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its own right, it is most prominent at high stehar masses 
and does not show a strong dependence on environment, 
so we do not consider it further here. 

4. DISCUSSION AND CONCLUSIONS 

Our results indicate a complex relationship between 
color, morphology, stellar mass, and group-centric dis- 
tance. Trends in color and morphology are distinct, and 
the use of a single property as proxy for galaxy type 
does not capture the whole picture. The most inter- 
esting trend seen is the shift in dominance among the 
low mass population from blue late disk galaxies in the 
outskirts to red bulge+disk types in group interiors. If 
part of one population was merely disappearing from the 
sample, either due to changes in mass or merging with 
other galaxies, then the other populations would all be 
expected to grow by an equal factor. The fact that the 
decline in one population is roughly balanced by the rise 
of a single other population suggests a transformation 
process. We now discuss these results in the context of 
past analyses, the physical implications of the present 
work, and future avenues to clarify the role of environ- 
ment in galaxy evolution. 

4.1. Connection to Previous Observations 

There is a long history of research into the covariance 
between the stellar masses, colors, morphologies, and en- 
vironments of galaxies and its evolution with time. Af- 
ter controlling for differences in stellar mass or luminos- 
ity, numerous studies at low redshift have found that 
color is more strongly correlated with environment than 
morphology (e.g., Kauffmann et al. 2004; Blanton et al. 
2005; Christlein & Zabludoff 2005; van den Bosch et al. 
2008; Bamford et al. 2009; Skibba et al. 2009; Weinmann 
et al. 2009). The implication of these studies is that 
the well-known correlation between morphology and en- 
vironment is secondary to correlations between morphol- 
ogy and stellar mass or color, with the latter properties 
more physically linked to environment. Still, some of 
these studies find residual correlations between morphol- 
ogy and environment after controlling for color and stel- 
lar mass or luminosity, particularly at low masses and 
among late to intermediate morphologies (Blanton et al. 
2005; Weinmann et al. 2009; Skibba et al. 2012). 

Moving from large low-redshift studies to moderate 
redshifts {z - 0.4), Balogh et al. (2009) and McGee et al. 
(2011) have shown that groups still host a higher frac- 
tion of red galaxies than the field, while Wilman et al. 
(2009) found an elevated SO fraction in groups with a 
mildly significant rise with group-centric distance. This 
radial trend in morphology is the opposite of that seen in 
Figure 3 but may be due to a smaller sample size or the 
use of luminosity- weighted centroids, which have been 
shown to be a poorer tracer of halo centers measured 
by weak lensing when compared with the most massive 
galaxy used in this catalog (George et al. 2012). 

Environmental trends at fixed stellar mass have been 
seen to weaken at higher redshift (e.g., Poggianti et al. 
2008; Tasca et al. 2009; Cucciati et al. 2010; lovino 
et al. 2010; Kovac et al. 2010), though recent analy- 
ses have shown clear differences between the fraction of 
star-forming galaxies across environments over a range 
of stellar masses at least up to z ~ 1 (Cooper et al. 
2010; Peng et al. 2010; George et al. 2011; Knobel et al. 



2012). Tasca et al. (2009) and Kovac et al. (2010) have 
compared the dependence of colors and morphologies on 
local density and in groups in the COSMOS field, find- 
ing a stronger effect on color than on morphology, sim- 
ilar to the low redshift results and suggesting a longer 
timescale for structural transformations than for quench- 
ing star-formation. Muzzin et al. (2012) and Presotto 
et al. (2012) have also shown that the quenched fraction 
among satellites depends significantly on group-centric 
distance at these redshifts. 

The results presented here extend this work to in- 
clude radial trends in colors and morphologies from 
z = 0.2 — 1. Our X-ray group catalog gives a large, clean, 
and fairly representative selection of 10^^ — lO"*^^ Mq ha- 
los (Finoguenov et al. 2010), and has been well-calibrated 
based on its weak lensing signal (Leauthaud et al. 2010; 
George et al. 2012). We detect significant trends in both 
color and morphology with group-centric distance. Some 
previous studies measured weak or insignificant gradients 
in morphology by using a simple dichotomy of spirals 
and ellipticals, and we can reproduce these results when 
using a coarse morphological binning (see discussion in 
Section 3.3.3). But in contrast to those results, we see 
clear morphological gradients at fixed stellar mass and 
color once the morphological classification considers dif- 
ferences in the bulge content of disk galaxies. 

The importance of this intermediate morphology be- 
tween pure disks and spheroidals has been noted before 
in the context of SO galaxies which have been known to 
dominate the evolution in the morphology-density rela- 
tion (Dressier et al. 1997; Postman et al. 2005; Smith 
et al. 2005; Boselh et al. 2006; Moran et al. 2007; Oesch 
et al. 2010; Lackner & Gunn 2013). Smith et al. (2005) 
also showed that morphological evolution occurs later 
in less dense regions than in the densest regions asso- 
ciated with clusters, with little evolution at intermedi- 
ate densities from z = 1 to 0.5 followed by an increase 
in early types at lower redshift. There is a significant 
overlap between the intermediate bulge+disk popula- 
tion in our analysis and visually classified SO galaxies 
although these bulge+disk galaxies are often classified 
slightly later along the Hubble sequence (Scarlata et al. 
2007). Our results in Figure 4 broadly corroborate those 
of Smith et al. (2005) and add that color and stellar 
mass are also important dimensions when studying these 
trends. 

4.2. Implications for Physical Mechanisms of Galaxy 
Transformation 

The radial gradients and redshift trends measured in 
Section 3 suggest a transformation among low mass 
satellites from blue late disks into red bulge+disks and 
spheroidals. The mechanism for this transformation 
must affect both color and morphology, or more phys- 
ically, the star formation rate and stellar kinematics. 
Some processes (see e.g. Boselli & Gavazzi 2006, for a 
review) halt star formation without significantly alter- 
ing stellar structure, such as gas removal via ram pres- 
sure stripping or weak tidal interactions, suppression of 
gas accretion within dense shock-heated environments, 
or quasar feedback. On the other hand, galaxy mergers 
and strong tidal interactions can affect both the distribu- 
tion of gas needed to form stars as well as the stellar mor- 
phology. A third possibility is that color and morphology 
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changes are physically coupled; bulge growth could sta- 
bilize a galaxy against disk fragmentation and suppress 
star formation (Martig et al. 2009), or gas loss may leave 
a disk unable to dissipate energy from tidal interactions 
or may drive instabilities leading to bulge growth. Yet 
another scenario is that a bulge only appears more promi- 
nent after a galaxy is quenched because the previously 
star- forming disk has faded. 

We can test these models by studying the morpholog- 
ical dependence on environment among quenched galax- 
ies. Mechanisms like gas stripping or disk fading that do 
not directly affect morphology should produce a higher 
fraction of quenched galaxies in dense environments, but 
among quenched galaxies the morphological distribution 
should be constant across environments. On the other 
hand, if quenched galaxies in dense environments have 
more bulge-dominated morphologies than quenched field 
galaxies, it would suggest that mergers or strong tidal 
interactions altering the structure of galaxies occur in 
addition to, or in conjunction with, the suppression of 
star formation. The results of this test are plotted in 
Figure 6, where we show the fractions of red galaxies 
that are early and late type as a function of environment. 
We have used a broad redshift range and combined the 
spheroidal and bulge+disk categories to reduce statis- 
tical errors for this smaller population of red galaxies. 
Figure 6 demonstrates that quenched satellites, particu- 
larly those in the inner regions of groups, are more likely 
to be bulge-dominated than their field counterparts. 

We interpret this to mean that some physical process 
in dense environments is driving bulge growth. Though 
some field studies (Bundy et al. 2010; Masters et al. 
2010) have already noted a tendency of passive disks to 
be more concentrated than star-forming disks, we find 
that the morphological evolution of quenched galaxies is 
even stronger in groups. Our result is consistent with 
previous arguments favoring bulge growth over disk fad- 
ing based on the higher typical luminosities of galaxies 
with intermediate morphologies compared to late types 
(Christlein & Zabludoff 2004; Burstein et al. 2005). The 
data presented here show that even at fixed mass there 
must be bulge enhancement that accompanies quenching 
in groups. 

The results of Figure 3 suggest that the timescale for 
the morphological transition from late disk to bulge+disk 
must be comparable to the timescale for quenching in or- 
der to turn blue late disks directly into red bulge+disk 
galaxies. For example, if quenching occurred via gas 
stripping or removal at a rate faster than any structural 
changes, we would expect to see blue late disks turn- 
ing into red late disks. Instead we observe a growth in 
red galaxies with earlier morphologies whenever blue late 
disks decline. The fraction of blue late disks plummets 
by more than a factor of two in ~ 2 Gyr among inner 
satellites (top left panel of Figure 4), while the fraction 
of red late disks is nearly constant and the fraction of red 
bulge+disks grows. Similarly, we do not detect a large 
population of blue bulge+disk galaxies which might be 
expected if bulge growth happened faster than quench- 
ing. These observations suggest that some morphological 
evolution occurs on a similar timescale as quenching. 

At the same time, the morphological transformation 
observed at low stellar masses has not fully turned many 
disks into spheroidals. The bulge+disk galaxies are still 
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Fig. 6. — Morphological fractions among red galaxies as a func- 
tion of distance from the group center. The small but significant 
excess in early-type morphologies among inner satellites relative 
to field galaxies suggests that mergers or tidal interactions cause 
more significant bulge growth among satellites. 

categorized as disks in the top level ZEST classification. 
This echoes earlier results finding a buildup of SOs in 
clusters but weaker evolution in the relative abundance of 
ellipticals (e.g., Dressier et al. 1997). Though the fraction 
of red spheroidals does not correlate significantly with 
group-centric distance (Figure 3), it does grow globally 
with time (Figure 4). This may indicate that different 
processes are responsible for producing bulge+disks and 
spheroidals. 

The discreteness of morphological classification makes 
detailed comparison of transformations in morphology 
and color somewhat difficult, but the structural evolution 
appears to contrast with that in color, where a clear bi- 
modality separates red and blue states and the relatively 
low intermediate population suggests a fast quenching 
timescale once it begins (e.g., Wetzel et al. 2012a, though 
see Balogh et al. 2011). One could test the hypothesis 
that galaxies pass through an intermediate state between 
blue late disk and red bulge+disk by measuring the time 
since quenching based on the SEDs of galaxies in these in- 
termediate states. Is there a difference in the mean color 
of blue bulge+disks and blue late disks, or between red 
bulge+disks and red late disks? Figure 5 does not show 
conclusive differences in the colors of these populations, 
and the satellite sample size is too small to constrain 
differences in the distribution of colors within the red 
and blue populations. More detailed analysis with spec- 
troscopic data could better constrain such evolutionary 
models. 

The group environment can play an important role 
in building up the well-studied color-morphology-density 
relations in nearby clusters. A group of mass 
at z = 1 should grow through mergers by an average of 
0.3 dex to z = (Fakhouri et al. 2010), so some of our 
high redshift groups will be the progenitors of massive 
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Fig. 7. — Illustrations of our main results and interpretation. Left panel shows the projected positions of satellites in an ensemble group 
with the same stellar mass and redshift range as Figure 2 with only the blue late disks (which dominate the outskirts) and red bulge+disks 
(which dominate the interior) displayed. The schematic diagram at right shows the effects of various physical mechanisms on color and 
morphology; the large arrow indicates the observed transformation from blue late disks to red bulge+disks, suggesting a combination of 
processes. 



clusters and others may accrete onto them. For compar- 
ison, a typical star- forming galaxy should grow by 0.2 
dex in stellar mass between each of our three redshift 
bins and by 0.7 dex from 2; = 1 to (Elbaz et al. 2011). 
Wetzel et al. (2012b) show that the dominant popula- 
tion of satellites in the most massive halos at z = were 
already satellites in groups when accreted. Our results 
demonstrate that both the color and morphology of these 
satellites are influenced in the group environment, show- 
ing the importance of "preprocessing" of galaxies prior 
to entering massive clusters. The bulge+disk population 
we see here is not quite as morphologically evolved as 
the growing SO population seen at low redshift and in 
massive clusters, but likely precedes it. 

We summarize our main results and interpretation in 
Figure 7. The left panel shows the positions of satel- 
lites for a stacked group and highlights the trend in Fig- 
ure 2 with blue late disks dominating the outskirts and 
red bulge+disk galaxies making up most of the inner 
satellite population. This combination of color and mor- 
phology transformations suggests some combination of 
gas removal and bulge growth, shown in the right panel. 
While the diagram is a simplification of a wide variety of 
model predictions, the primary point is that both color 
and morphology are affected in group satellites, so phys- 
ical mechanisms that explain environmental correlations 
should refiect this. 

4.3. Future Prospects 

While these measurements add insight to the transfor- 
mation of galaxy properties, much work remains to be 
done to constrain the variety of physical mechanisms that 
may be responsible. A cleaner bulge-disk decomposition 
would be useful to track the growth of bulges and the 
importance of disk fading (e.g., Lackner & Gunn 2012, 
2013). Incorporating the physical sizes of each compo- 



nent would also help to link these transformations with 
the significant growth seen in early types since z ^ 2 
(e.g., Bruce et al. 2012). With deep high-resolution imag- 
ing from HST, the CANDELS survey (Grogin et al. 2011; 
Koekemoer et al. 2011) is pushing these studies to higher 
redshift and should also allow bulge-disk decomposition 
to be studied at multiple wavelengths, identifying when 
and where star formation is happening or stopping. Al- 
ternatively, visual classifications from an ongoing Galaxy 
Zoo project in the CANDELS fields will track the evo- 
lution of spiral features and bars which are challenging 
for automated techniques. Impending wide field imaging 
surveys will add greatly to the statistics of these analyses. 

Different measures of environment may also help dis- 
entangle the relevant physical mechanisms. While we 
have argued in favor of halo-based indicators (see George 
et al. 2011 for a discussion), the local galaxy density 
within a halo^ in addition to group-centric distance, may 
shed light on galaxy-galaxy interactions and infalling 
substructure (e.g., Blanton & Berlind 2007; Cibinel et al. 
2012; Woo et al. 2012). However, this indicator is a noisy 
quantity affected by shot noise, redshift errors, and pecu- 
liar velocities, likely requiring deep and complete spectro- 
scopic data for clean results. Larger surveys will also en- 
able a wider range of halo masses to be probed at higher 
redshift. Comparison of different halo mass proxies can 
provide a test of environmental mechanisms, for example 
if X-ray bright groups are more efficient at ram pressure 
stripping satellites than groups with less hot gas. 

These studies can also be extended to the interplay of 
environmental mechanisms with non-stellar components 
of galaxies, namely the gas content and central black 
holes in galaxies. Current and planned radio arrays will 
extend the study of neutral and molecular gas beyond 
the local Universe allowing a clearer picture of how star- 
formation is fed and quenched. And the growth of bulges 
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seen in this paper should be connected with accretion 
onto the central black hole in order to explain the tight 
correlation seen locally between these components. 

Finally, while we have simply presented an empirical 
description of the data here, we can extend this work by 
modeling the evolution of color and morphology with en- 
vironment. Peng et al. (2010) and Wetzel et al. (2012b) 
present simple empirical descriptions of the fraction of 
quenched galaxies as a function of stellar mass and envi- 
ronment. Incorporating morphologies into this frame- 
work will further illuminate the physical processes at 
work that build up the environmental correlations we 
observe. 
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